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SUMMARY 


A reliable and efficient iterative method has been developed 
for designing wing section contours corresponding to a prescribed 
subcritical pressure distribution. The desiyn process is ini- 
tialized by using MCAERO (MCAIR 3-D Subsonic Potential Flow Analy- 
sis Code) to analyze a baseline configuration. A second program 
DMCAERO is then used to calculate a matrix containing the partial 
derivative of potential at each control point with respect to 
each unknown geometry parameter by applying a first-order expan- 
sion to the baseline equations in MCAERO. This matrix is calcu- 
lated only once but is used in each iteration cycle to calculate 
the geometry perturbation and to analyze the perturbed geometry. 
The potential on the new geometry is calculated by linear extra- 
polation from the baseline solution. This extrapolated potential 
is converted to velocity by numerical differentiation, and velo- 
city is converted to pressure using Bernoulli's equation. There 
is an interactive graphics option which allows the user to 
graphically display the results of the design process and to 
interactively change either the geometry or the prescribed pres- 
sure distribution. Not only is this design procedure accurate 
for large perturbations, but the cost of each iteration cycle is 
more than two orders of magnitude less than a conventional 
analysis solution. Examples of the design process are presented 
to demonstrate that the method is accurate, numerically stable, 
and converges in only three to five iterations. 
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INTRODUCTION 


The surface singularity panel method has been widely 
accepted as an excellent means of determining the subcritical, 
potential flow about complex aircraft configurations (References 
1-9). The better formulated methods will consistently predict 
accurate wing pressure distributions, even in regions with strong 
fuselage-nacelle-store interfernce. A typical example of this 
accuracy obtained by MCAERO (MCAIR 3-D Subsonic Potential Flow 
Analysis Code) is shown in Figure 1 for an F-18 configuration. 
The power of surface panel methods has been recognized by several 
investigators who have developed iterative inverse techniques for 
designing wing design section geometry corresponding to 
prescribed pressure distributions (References 10-13). However, 
each of these investigators have encountered and not overcome at 
least two of the shortcomings listed in Figure 2. 



Angle-of*Attack - deg Pitching Moment Coefficient, C M 
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Figure 1. MCAERO Prediction Accuracy for Clean and Takeoff Configurations 
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Symptoms 

Cause 

Remedy 

o Calculated Pressures Overly 
Sensitive to Contour Smoothness 

o Source Singularities 

o Combined Source - Doublets 
(Green’s Identity) 

o Calculations Do Not Converge 

o Inaccurate Leading Edge 
Design Contours 

o Dirichlet Boundary Conditions 

© Transpiration Principal 
(Equivalent Blowing) 

© First-Order Mathematical 
Expansion (3C p /3z) 

« Unrealistically Wavy 
Design Contours 

o One Unknown Coordinate per 
Prescribed Pressure 

o More Pressures Than Unknowns 
o Least Squares 

o High Cost per Iteration 

o Complete Panel Method Analysis 
During Iteration Cycle 

o Perturbation Analysis Method 
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Figure 2. Common Barriers to a Successful Wing Design Method 


A reliable and extremely efficient method for solving the 
wing-on-fuselage design problem depicted in Figure 3 was 
developed in 1982 (Reference 14). The development of this method 
by McDonnell Aircraft Company (MCAIR) was supported under 
contract to NASA Langley Research Center (LaRC) . The efficiency 
is the result of employing the pertubation analysis method of 
References 8 and 9 in each iteration cycle of the design 
solution. The design approach is based upon the execution of 
three complementary but independent computer programs. For a 
given baseline panelled configuration, the first program 
generates a conventional panel method analysis solution, 
including surface singularity strengths and pressure 
distribution. The solution file from the first program provides 
the input to the second program, which calculates a matrix of 
partial derivatives of surface potential with respect to 
arbitrary geometry perturbations . The first and second program 
are executed one time only in order to provide a permanent input 
file to the third program, which calculates the winy geometry 
corresponding to a prescribed pressure distribution. 

As part of the MCAIR 1982 and 1983 Independent Research and 
Development Project, user-oriented production versions of the 
above three computer programs were developed. These programs are 
respectively designated "MCAERO" , "DMCAERO" , and "DESIGN”. A 
fourth production program that applies the perturbation analysis 
method of References 8 and 9 is designated "PAM" . 

As part of the present contract with LaRC (NAS1-17176 ) , an 
interactive graphics option for the design method was estab- 
lished. This report summarizes the capabilities of the user- 
oriented production codes and the interactive graphics design 
option . 
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Given: o Fuselage geometry 

o Baseline wing 
o Angle-of-attack 

o Prescribed pressure distribution (C p .) 
Calculate: Wing section geometry 



Use of trade names of manufacturers in this report does not 
constitute an official endorsement of such products or manufac- 
turers, either expressed or implied, by the National Aeronautics 
and Space Administration. 
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2. PERTURBATION ANALYSIS METHOD 


The standard solution approach to prescribed pressure design 
problems is to divide each iteration cycle into an analysis, 
pressure calculation step and an inverse, geometry correction 
step. In the two-dimensional method of references 7 and 15, an 
entire panel method solution is calculated during each analysis 
step. Furthermore, a new geometry-velocity perturbation matrix 
is calculated for each inverse step. In spite of the fact that 
the number of computations in each iteration cycle is a cubic 
function of the number of panels, the total computing cost is 
relatively small. The reason is that with the better formulated 
methods typical two-dimensional problems require fewer than one 
hundred panels. 

However, the number of panels required for wing-fuselage 
configurations is almost an order of magnitude greater, and the 
approach of the two-dimensional design procedure would be 
extremely expensive. During the initial development of the pilot 
code version of the wing design method, it became evident that 
the cost of a practical wing-on-fuselage design procedure would 
be prohibitive using the then existing panel method technology. 
Therefore, in order for the design procedure to become practical, 
the existing panel method technology was expanded by a new cost 
savings feature, i.e., the Perturbation Analysis Method (refer- 
ences 8 and 9). The method is an extremely efficient tool for 
analyzing the pressure distribution corresponding to a series of 
arbitrary, small perturbations to a baseline wing— fuselage geo- 
metry (Figure 4). The following features of the perturbation 
analysis method make it especially practical for application to 
an iterative wing section design method. 

(1) The computational expense for analyzing each successive 
geometry perturbation is almost two orders of magnitude 
less than that of a conventional panel method, 

(2) The pressure distribution prediction accuracy is compe- 
titive with conventional surface panel methods for very 
large perturbations to wing section geometry. 

(3) A pre-calculated matrix of partial derivatives for the 
paneled baseline configuration is available. Each ele- 
ment of the matrix is the rate of change of potential 
at a boundary condition control point with respect to a 
geometry parameter perturbation. For design applica- 
tions, the geometry-potential perturbation matrix can 
be efficiently converted to a geometry— pressure pertur- 
bation matrix. 

The purpose of this section is to give a brief review of the 
fundamentals and power of the Perturbation Analysis Method (a 
more complete description is available in references 8 and 9). 
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Objective 

o Subsonic inviscid analysis of multiple geometry 
perturbations at small additional cost 

Approach 

o Precalculated baseline matrix of potential 



GP23-0838-5 

Figure 4. Perturbation Analysis Method 


2.1 MATHEMATICAL APPROACH - The perturbation analysis approach 
requires an initial baseline calculation of a conventional panel 
method solution for an arbitrary baseline configuration. Subse- 
quently, a matrix consisting of the analytically derived partial 
derivatives of velocity potential with respect to geometry coordi- 
nates is calculated. The baseline solution and derivative matrix 
are calculated one time only and then stored for repetitive use. 
For each geometry perturbation, the solution surface distribution 
of velocity potential is constructed by multiplying the deriva- 
tive matrix by a new right-hand-side. This procedure bypasses 
the two computationally expensive steps of a conventional panel 
method: calculating the influence coefficients and solving a 
large system of linear algebraic equations. 

Although the perturbation analysis method is appropriate for 
predicting the effect of arbitrary small changes to wing planform 
and fuselage geometry, the real power of the method is the accu- 
racy with which large perturbations to wing thickness, camber, 
and twist can be analyzed. Examples of the accuracy and effi- 
ciency of this method can be found in references 8, 9, and 14, 
and therefore will not be repeated here. 
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The first three computer codes listed in Figure 5 are 
required in order to use the production version of the Perturba- 
tion Analysis Method. The first program, MCAERO , is the MCAIR 
conventional panel method for analyzing the baseline configura- 
tion. The second program DMCAERO, employs a differential 
mathematical formulation and an output file from MCAERO to 
calculate the matrix of partial derivatives of the perturbation 
potential. For each baseline configuration, DMCAERO creates an 
input file for the third program - Perturbation Analysis Method 
(PAM). The pilot code versions of both MCAERO and DMCAERO were 
developed basically to validate the procedure. While both worked 
successfully, each suffer from a similar problem, tedious input. 
Therefore a prime concern in the development of the production 
version of the codes was to simplify the input while improving 
efficiency and accuracy. Figure 6 shows the success obtained in 
simplifying the input for both MCAERO and DMCAERO. 



MCAIR 3-D 
Subsonic Potential 
Flow Program 
(MCAERO) 

Geometry Influence 
Coefficient 
Program 
(DMCAERO) 

MCAIR 3-D 
Perturbation 
Analysis 
Program 
(PAM) 

MCAIR 3-D 
Wing Design 
Program 
(DESIGN) 

Input 

• Baseline 
Geometry 
(x, y, z)j 

© Baseline 
Geometry 
(X, y, z)j 

© Baseline 
Potential 
Distribution 

o Baseline 
Properties 
(x. y, z),, <t> j 

© Derivative Matrix 
(d<t>/dx:, 90/9^, 
dfydzj 1 

© Geometry 
Perturbation 
(Ax, Ay, Az)j 

© Baseline 
Properties 
(X, y. z)j 0 i 

• Derivative 
Matrix 
(90/a Zj) 

• Prescribed 
Pressures 

Approach 

• Conventional 
Panel Method 

o First-Order 
Expansion to 
Panel Method 
Formulation 

• Linear 
Extrapolation 
A0,= 

W/dXj) Ax, 
i 

+ (30/9^) Ayj 
+00/3 z,) Azj] 

• Linear 
Extrapolation 
and Least 
Squares Solution 

Output 

• Baseline 
Potential 
Distribution 0 ( 

© Baseline 
Aerodynamic 
Properties 

• Derivative 
Matrix 
(90/9x,, 
H/dy-y 

• Aerodynamic 
Properties 
of Perturbed 
Geometry 

• Designed 
Wing Geometry 
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Figure 5. The Computer Programs for the 3-D Wing Design Method 
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Number 

of 

Inputs 
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Figure 6. Number of Inputs for a Typical Case 


While MCAERO and DMCAERO can be relatively expensive to run, 
once the input file for PAM has been generated, the first two 
programs are no longer required. PAM can be executed repeatedly 
at low computing cost for the analyses of a series of perturba- 
tions to the panel corner coordinates (x,y,z)j of the complete 
aircraft configuration. The method used to analyze each perturba- 
tion is the same as the conventional panel method calculation 
with two significant exceptions. First, no influence 
coefficients are calculated; second, no large system of linear 
algebraic boundary condition equations is solved. Instead, the 
perturbation potential at each control point is calculated by 
linear extrapolation. The conversion of potential to surface 
velocity is based on numerical differentiation. 

Usually the second program is executed twice for each base- 
line configuration, once at 0° incidence and once at 90° inci- 
dence. By employing the principle of linear superposition, the 
perturbation analysis is automatically performed at any interme- 
diate angle of attack (Figure 7). 
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Figure 7. Linear Superposition of 0° and 90° Solutions 
3-D Wing Design Method 


Figure 8 shows the relative cost of each of the codes listed 
in Figure 5 for a typical wing only and wing-on-fuselage case . 
As can be seen once the "expensive" codes (MCAERO and DMCAERO ) 
are run, the computing costs become insignificant. 

W = Wing only typical case (200 
WF = Wing-fuselage typical case 

300 10 


8 

200 

6 

Relative 

Computing 

Cost 

4 

100 


2 


0 0 



w 









i 




i 



1 


w 

I 


I 


I— MCAERO— 


-DMCAERO-] 


panels) 

(700 panels) 


WF 



Figure 8. Relative Computing Costs for the Four Codes Used in the 

Wing Design Method 
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3. 3-D WING DESIGN METHOD 


For a given panelled wing or wing— fuselage configuration, 
assume that prescribed pressure coefficients for a fixed arbi- 
trary angle of attack are assigned to the panel centers of the 
wing* The objective of the design method is to determine the 
change in wing section geometry that most nearly corresponds to 
the prescribed pressure distribution (Figure 3). 

As was shown previously, the Perturbation Analysis Method 
can accurately calculate the pressure distribution corresponding 
to large changes in wing section geometry. With that as a basis, 
it was reasoned that an inverse formulation to the Perturbation 
Analysis Method could accurately design large changes in wing 
section geometry. It was this approach that formed the basis of 
the wing design method. The design method must employ an itera- 
tive scheme because pressure is a nonlinear function with respect 
to geometric changes. However, recognition of the fact that the 
perturbation potential is a nearly linear function eliminated the 
need to perform extensive computations in each iteration cycle. 

A schematic of the wing design method is presented in Figure 
^ * The method can be used to solve very general aerodynamic 
design problems. For example, the prescribed aerodynamic quan- 
tity at a panel center can be either a velocity component or 
pressure coefficient. Arbitrary geometry parameters such as wing 
chord or fuselage shape can be selected for design. Furthermore, 
design constraints such as fixed camber or thickness can be 
imposed. Most design problems, however, are of the type illus- 
trated in Figure 3. This type, designated the "standard wing 
design problem", is defined in detail below. The remainder of 
this section presents the mathematical formulation for the wing 
design method and example design solutions . 

3.1 STANDARD WING DESIGN PROBLEM - As illustrated in Figure 10, 
the region of panels subject to design is identified by corner 
points in the range (^A'Da) _< (i»j) £ (iB'jfi)' where the limits 
( 1 a/ iB» jB) are selected by the user. If jA and je are points 
on the upper and lower surface trailing edge respectively, then 
the geometry of the complete wing section at each span station i 
will be designed. One restriction of the wing panelling is that 
it be basically trapezoidal. At the center of each panel in the 
design region, the desired pressure coefficient is prescribed by 
the user. 
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Given: Baseline Configuration Panelling 

_ 84); 

Precalculated Arrays: C n , and . 

' az k 


5 


User Input: Prescribed Pressure Distribution 

K ^ 

SI 

r 

Convert Matrix [dfald z k ] to [3C p ./3z k ] 

4 



Solve Linear, Alge 

A °p, e E 

k 

Oraic Equations for Az k 

5C Pil A 
• Az k 

3z k J 


[I 


Update Geometry: z k = z k + Az k 


J 


Update Pressure: Perturbation Analysis Method 


X 

-i 

^0 

Test for Convergence 

1 1 1 


" ST 

Solution 


G P 33 0 692 7 

Figure 9. Schematic of 3-D Wing Design Method 
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The unknowns are Az at the panel corner points {(iA»jA + l) 
_< (i»j) £ * However, not all of the unknowns are per- 

mitted to be independent. As illustrated in Figure 11, less than 
one-half of the unknowns are independent. The remaining unknowns 
are generated by interpolation through the independent unknowns . 
On the span stations i = i^, i&+2, iA+4, ••., iB# each dependent 

unknown Az(i,j) is established by least squares quadratic interpo- 
lation through the path of points (j-3, j-1, j, j+1, j+3 ) on the 

baseline configuration. For the remaining span stations (i = 
iA+1/ iA +3 * . .., iB -1 )/ each unknown Az(i,j) is established by 

straight line interpolation through Az(i-l,j) and Az(i+l,j). For 
this type of interpolation to be accurate, it is implicitly 
assumed that the three points (i-l,j), (i,j), and (i+l,j) lie on 

nearly the same per cent chord line. Typical wing panelling is 
consistent with this assumption. 
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“Even” Span Stations (i A , i A + 2,..., i B — 2, ig) 

For each “even” j (j A + 2. J A + 4 JB - 4 , IB " 2 > 



“Odd” Span Stations (j A + 1. i A + 3 i B — 3, i B — 1) 

Az (i, j) is generated by linear interpolation between 
span stations i - 1 and i + 1 


GP33 0151-14 

Figure 11. Interpolation Scheme 
for Dependent Unknowns 
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The reason for limiting the number of independent unknowns a 
priori is to prevent numerically unstable design calculations . 
Figures 12a and 12b illustrate two types of design instabilities 
that could occur if every panel corner in the design region were 
allowed to be an independent unknown. 


a) Wavy wall instability 

Pressure at panel center does not control 
streamwise slope oscillations 



b) Four-corner instability 

Average plane of panel does not control 
corner point oscillations 



GP330151-13 

Figure 12. Typical Wing Design Solution Instabilities 


Consistent with the nomenclature of the Perturbation Analy- 
sis Method (reference 8), each independent unknown perturbation 
is assigned an index kg. The value of Az for perturbation number 
kg is designated S^g. A schematic of the independent unknowns is 
presented in Figure 0 13. The objective of the wing design method 
is to calculate the values Sj^g (1 £ ks £ NKS ) such that the 
resulting pressure distribution most nearly corresponds to the 
prescribed pressure distribution. 
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Figure 13. Ordering of Independent Unknowns S ks 


3.2 MATHEMATICAL FORMULATION FOR WING DESIGN METHOD - For any 
panelled baseline configuration, application of the wing design 
method requires that the arrays ^ and 3(^/3 S^ g have been calcu- 
lated a priori. <t> i is the perturbation potential at the ith con- 
trol point and is the rate of change of <f> -j_ with respect 
to independent geometry perturbation number kg. As described 
previously MCAERO and DMCAERO will automatically calculate the 
required arrays and store them on a computer disk file. The wing 
design method can then be used to calculate the geometry perturba- 
tion that most nearly matches prescribed aerodynamic properties 
within the limitations of a minimal least square error. 

At any panel center selected by the user, one or more proper- 
ties can be prescribed. The property can be either pressure 
coefficient (Cp) or velocity component in an arbitrary, specified 
direction. The prescribed value of an aerodynamic quantity at a 
panel center is designated Qp^, where there is one index i for 

each prescribed value ( 1 < i _< Njjes)* 
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Arbitrary geometric constraints can be imposed upon the inde- 
pendent unknowns (S^g). Each geometric constraint is expressed 
as a linear equation 

n d 

Z Cc iD ' ^sUd ) 3 = CRHS U) 

i D =i 


Where Np, ci D , and CRHS are arbitrary values specified by the 
user. Depending upon the values specified, the constraint 
equation can be used to fix the cross-sectional area of a wing 
section, fix the thickness at one point, and so forth. 


The aerodynamic design problem can now be expressed in mathe- 
matical form. The objective is to calculate the array of indepen- 
dent geometric unknowns S^g (1 _< kg <_ NKS) that will minimize the 
following function E. 


E 



AREAj_ • 



N CON 

♦ E 

3=1 


WCON 


2 . 


n D 

c ip * s ks(ip) 

i D =! 


2 

CRHS 

J j 


( 2 ) 


where (1) 


each i (1 _< i _< N DE g) corresponds to one prescribed 
aerodynamic quantity at one panel center. 


(2) each j (1 <_ j £ n C0n) corresponds to a geometric con- 
straint equation. 


(3) the weights WDES-j_ and WCONj are specified by the user 
(typically, WDES^ = 1.00 and WCONj >>1.00), 

(4) AREAi is the area on the baseline configuration of the 
panel corresponding to prescribed quantity Qp^» and 

(5) (Qi - Qp . ) is the difference between the calculated 
and prescribed values of an aerodynamic property. 


Iteration is used to solve for the array of unknowns S^g 
corresponding to a minimum value of E from equation 2. As shown 
in Figure 9, each iteration cycle is divided into an inverse step 
in which the geometry perturbation is calculated and a direct 
step in which the perturbed or updated geometry is analyzed. In 
order to reduce costs of the design method, the pressure distribu- 
tion for the complete configuration is analyzed only for the 
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zeroth (baseline) and final iteration cycle. For the intermedi- 
ate cycles, the pressure distribution for only the design region 
is calculated using a procedure depicted in Figure 14 . As 
illustrated, the spanwise and chordwise components of the 
perturbation velocity (designated and qs) are calculated using 
a one-dimensional interpolation through (Uj_, V*3 and y 5 ) and (U2» 

V*3, and U4), respectively. Once qa. and qe are known, the total 
velocity vector for the panel is determined as follows 




{q A ; q B } is {Spanwise; Chordwise} Component of 
Perturbation Velocity and is Calculated 
by One-Dimensional Interpolation Through 

{(Mil M31 M5); (M2 ’ M3> M 4 )}- GP33-0692-4 


Figure 14. Velocity Calculation on Design Panels 
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( 3 ) 


$ = {\? + ae. T + q_ e_ + q e } 

oo N B c 

where 

is the panel source strength 
e^j is the panel normal vector 
q c is defined as 

„ _ 9A ~ <lB e B * e A 

4c - ~ — 

e c * e A 

Sc 1S e B x e N 

In the inverse step of each iteration cycle, the matrix of 
derivatives is calculated. Then the change in Qi induced by a 
small perturbation to the array S^g can be expressed as 

NKS 

V 9U i 

dQi = ^ -as]- ' ds *s (4) 

k s =i k s 

By incorporating equation (4) in equation (2) and minimizing E 
with respect to dS k g, a system of linear, algebraic equations is 
established. The solution by standard matrix algebra provides 
the values dS^g. The updated geometry is then analyzed by the 
Perturbation Analysis Method. 

3Ui 

One might expect that recalculation of the matrix t-s— 3 

dSk s 

during each iteration cycle would require substantial computing 
expense. However, the following approach has proved to be both 
very efficient and accurate. 

For each panel in the design region, an aerodynamic quantity 
Qi, which can be either pressure coefficient or a velocity 
component at the panel center, is calculated. If Qi is a 
pressure coefficient then 


Qi - c p 2/V RE p [Vco x V x + Vco y V y + Voo z V Z J 

- C6 V x + Vy + V z 3 /V ref 

3 2 2 

Jg^— = -2/Vref [Vc» x + 6 V x , V OTy + Vy, V OOZ + V z 3 



( 5 ) 


(b) 
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where 


8 2 = 1 - M 2 

H 00 

v x» v y» V 2 are the perturbation velocity components. 

If Qi is a velocity component in an arbitrary fixed direc- 
tion (with direction cosines COSX, COSY, COSZ), then 

Qi = V • (COSX, COSY, COSZ) (7) 

9Qi a 

3Sj^ = (C0SX ' C0SY ' C0SZ) * CV X' V (8) 

As shown by equations (6) and (8), the problem of calculat- 
ing the derivatives of Qi essentially reduces to calculating the 
velocity derivatives 

9V X 3 V y 9V Z ^ 

9S k S ' 9S ks ' 9S ks ) 

At the center of the panel the velocity vector V can be 
represented as 

->■ ->■ 

V = V + V(j> (9) 

The gradient of the perturbation potential V<f> can be calcu- 
lated by numerical differentiation through 21 neighboring control 
points j. Mathematically this is 

21 

V<f> = 2 <$> j * ( a j e x + b j e y + c j e z ) (10) 

j=l 

where the scalars a j , bj and cj are functions of the geometry at 
the panel corner points . 

By substituting equation (10) into equation (9) and differen- 
tiating, the desired velocity derivatives can be expressed as 


3 S kg ( V x* v y/V z ) - 


+ 


21 


^ 9S k * (a j e x + b j e y + c j e z) 




( 11 ) 


4* j ’ 3s kg (a j e x + b j e y + c j e z) 


19 



The only term in equation (11) that requires substantial expense 

to compute is Fortunately, it is also the only term that 

° 

is nearly independent of perturbations to wing thickness, camber, 
and twist. Therefore, the precalculated baseline matrix C-|o~3 

dSfcg 

that is available on computer disk file from DMCAERO can be used 
in equation (11) during every iteration cycle. Substitution of 
equation (11) into equation (6) or (8) yields the desired value, 
9Qi 


A significant feature of the preceding approach is that the 

accuracy of the calculated matrix .. 9Qi ... . . 

J is competitive with an 


exact first order expansion during each iteration cycle. How- 
ever, much less computing effort is required. 


In fact, the number of computations required for a complete 
iteration cycle is relatively small. The reason is apparent upon 
consideration of each calculation step in Figure 9. For example, 
consider the system of linear, algebraic equations to be solved 
for the perturbations to the independent unknowns. Typically, 
fewer than one hundred unknowns are sufficient for wing design, 
compared to several hundred for a conventional panel method 
analysis of a wing-fuselage. Also, consider the last step of 
each iteration cycle - analyzing the updated geometry. The 
extremely efficient Perturbation Analysis Method is used for that 
calculation. This coupled with the fact that for 'the interme- 
diate cycles only the design region is analyzed results in a very 
efficient method for designing wing section geometry. 

The wing design method has been automated and is fully opera- 
tional on the McDonnell Douglas CYBER 176. The production 
computer program is designated '’DESIGN". A demonstration of the 
accuracy, efficiency, and numerical stability of the method is 
presented in the next section. 

3.3 EXAMPLE DESIGN SOLUTIONS - Two example solutions by Program 
DESIGN are presented below. The calculations were performed on 
the McDonnell Douglas CDC CYBER 176. 

In the first example, the objective is to design a 
fighter-type wing section geometry by starting from the baseline 
NACA-0012 wing panelling of Figure 15. The fighter wing pressure 
distribution at 0° angle of attack was prescribed at the center 
of each of the 208 panels. The converged solution after 3 
iterations is presented in Figure 15. The relative cost of the 
designed solution is shown in Figure 8 where the design solution 
costs approximately 1/4 that of the conventional MCAERO solution. 
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Baseline wing section (NACA 0012) 

Target geometry (fighter) 

O Converged solution (3 iterations) 


a) Geometry 



b) Pressure Distribution 


QP23-0S3S-S 


Figure 15. Example Wing Design Solution 





In the second example, the objective is to design a wing on 
a fuselage-nacelle-tail configuration. The baseline is shown in 
Figure 16. For this case, the prescribed pressure was basically 
a linear pressure distribution which removed the leading edge 
pressure spike. The resulting pressure and geometry are shown in 
Figure 17. Although the solution converged in two iterations, a 
third iteration was calculated for verification. The computing 
cost of the design solution was only about 1/25 that of the 
conventional solution (see Figure 8). Based upon these and other 
test cases, it was concluded that the 3-D Wing Design Method is 
indeed an accurate, efficient and cost effective method of 
designing wing section geometries corresponding to prescribed 
pressure distributions. The next section describes the 

interactive graphics module. 


Full 3-D Panelling 

• 596 Panels Per Side of Symmetry Plane 

• Flow-Throuh Nacelles 
® Panelled LEX 

• Panelled Horizontal and Vertical Tails 



GP330692-3 


Figure 16. Isometric of Panelled F-18 
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4. INTERACTIVE GRAPHICS OPTION 


The interactive graphics option is a stand-alone self- 
documented code which takes an output from the 3-D wing design 
code, plots user selected pressure and geometries, and if 
desired, creates a new input file, interactively, for the design 
code. This graphics module has 4 options built into it. The 

first allows the user to plot any or all the geometry segments in 
a three view plot (Figure 18). The second option allows the user 
to plot a wing section geometry and the corresponding pressure 
distribution (Figure 19). 





GP330692-13 


Figure 18. Three View Geometry 

Interactive Graphics 
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x/c 

Plot Label 
Y = 151. 15000 

Display Additional Geometry and C p Distr (Y/N)? gp33-069214 


Figure 19. Wing Section and Pressure Distribution 
Interactive Graphics 


Options 3 and 4 are used when the user wants to modify 
either the initial wing geometry or prescribed pressure distribu- 
tion for another design case. Option 3 allows the user to modify 
the initial guess of the starting geometry. The code displays 
the selected wing section, both initial baseline and the result- 
ing designed wing (Figure 20). If the user requests geometric 
changes, then the point numbers are placed on the points that the 
user may change (Figure 21), the screen is then cleared, and the 
point number and value are displayed (Figure 22). The user may 
then select new values. When done, the code displays the new 
geometry (Figure 23). 
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Change Next Iteration Geometry? (Enter Y or N) 


— Baseline 


^ Calculated 



X 


Plot Label 
Wing Section = ? 

Y = 181.40000 

GP33-0692-1 5 

Figure 20. Wing Section Geometry 

Interactive Graphics 


Change Next Iteration Geometry? (Enter Y or N) 
?Y 

Type “C” to Continue 
?C 


— — Baseline 


^ Calculated 



X 


Plot Label 
Wing Sections? 

Y = 181.40000 

GP330692-16 

Figure 21. Wing Section Geometry 
Interactive Graphics 
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IINT I 

1 VALUE 

1 

86 • 34609 

2 

86.66900 

3 

86.15800 

4 

84.95500 

5 

83.74300 

6 

82.97300 

7 

82.38900 

8 

81.96400 

9 

81.98400 

10 

83.22100 

11 

85.79700 

12 

86.30700 


ENTER POINT NO.,NEU 2 VALUE 
ENTER 0,0 TO END 

71.87 

72.87 
^11 ,87 

712.87 
70 0 

DISPLAY GRID? 

GP330692-17 

Figure 22. Wing Geometric Points 

Interactive Graphics 


Accept Changes to This Section (Enter Y or N) 
?Y 



GP33-0692-18 


Figure 23. Updated Wing Section Geometry 

Interactive Graphics 
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By using Option 4, the user may change the prescribed pres- 
sure distribution. Initially the code displays both the calcu- 
lated and prescribed pressure at the requested station (Figure 
24). If the user wants to change any pressures, the pressure 
points are numbered (Figure 25), screen is cleared, and the point 
numbers and pressures are displayed (Figure 26). The user may 
then select the desired points and their new values. When done, 
the code displays the new prescribed pressure as shown in Figure 
27. Once the user is satisfied with the new pressure and geo- 
metry, the code generates an input file for the wing design code. 
This iterative procedure can be repeated again until a final 
desired solution is obtained. 


Change Prescribed Pressures? (Enter Y or N) 
? 

— Calculated C p 


O Prescribed C p 



Figure 24. Wing Section Pressure Distribution 
Interactive Graphics 
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Change Prescribed Pressures? (Enter Y or N) 
?Y 

Type “C” to Continue 
? 


Calculated C p 


0 Prescribed C p 



Plot Label 
Y= 151.15000 


GP33-069220 

Figure 25. Wing Section Pressure Distribution 

Interactive Graphics 


POINT S 

PRESSURE 

1 

0*00000 

S 

-*20000 

3 

-*23000 

4 

-.27000 

5 

-.33000 

6 

-.40000 

7 

-.59000 

8 

-.75000 

9 

-1.25000 

10 

-1.60000 

11 

-1.70000 

12 

-1.60000 

13 

-.10000 

14 

.75000 

15 

.64000 

16 

.50000 

17 

.25000 

18 

-.08000 

19 

-.18000 

20 

-.12000 

21 

.13000 

22 

.33000 

23 

.42000 

24 

.36000 

ENTER POINT N0.,NEU CP UAUUE 
ENTER 0,0 TO END 
?6, .45 
?7, .65 
?8, .90 
?0,O 

Figure 26. Wing Pressures 

GP 33-0692-21 


Interactive Graphics 
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Accept New C p Distribution for This Chord-Wise Strip? (Enter Y or N) 
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Figure 27. Modified Wing Pressure Distribution 

Interactive Graphics 
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5. CONCLUSIONS 


The perturbation analysis and wing design methods are simi- 
lar to ^classical thin wing theory in the sense that small distur- 
bance . linearized assumptions are employed. This mathematical 
simplification generates extensive computational savings for aero- 
dynamic problems involving successive iteration, such as design. 

On the other hand, the restrictions of classical thin wing 
theory have been eliminated. Compared to an exact potential flow 
solution, the present approach is quite accurate for thick wings, 
large leading edge radius or camber, and high angle of attack. 
The success of the wing design method is attributed to the inclu- 
sion of all significant first-order geometry-pressure perturba- 
tion terms in each iteration cycle. This leads to rapid solution 
convergence, in spite of the fact that the entire distribution of 
surface potential is constructed by simple linear extrapolation. 
With the inclusion of the interactive graphics module, the wing 
design method has become a unique, powerful, highly efficient 
method for designing wing geometries which correspond to 
prescribed pressure distributions. 
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